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We study the action and the energy densities of the confining string in the indirect Z{2) projection of SU(2) 
lattice gauge theory. We find that the width of the confining string is proportional to the logarithm of the distance 
between the quark and antiquark. In 2q2q system we observe the effect of the reconstruction of the flux tube 
when we change the distance between two qq pairs. 



1. INTRODUCTION 

The structure of the confining string in terms 
of the Abelian or monopole operators after the 
Abelian projection has been studied in lattice 
SU(2) and SU(3) gluodynamics and in lattice 
QCD with two quark flavors 1 2 . The distri- 
bution of the monopole currents and the electric 
field near the confining string is similar to the 
distribution of the currents of the Cooper pairs 
and magnetic field near the Abrikosov string [3]. 
Thus the quantum confining string is similar to 
the classical solution of the equations of motion 
of the Abelian Higgs model. For the center- 
vortex model of confinement 0] there is no corre- 
sponding classical model, confinement is due to 
(quasi)randomly distributed magnetic fluxes in 
the vacuum. On the other hand the numerical 
investigation of the confining string in terms of P— 
vortices is even simpler than that in the Abelian 
projection, since the statistical fluctuations are 
smaller In Section 2 we give the main def- 
initions. In Section 3 we present the results of 
the calculations of the electric and magnetic fields 
near the confining string. In Section 4 we discuss 
the confining strings in the 2q2q system. 

All results are obtained on 24 4 lattice at = 
2 AO (20 statistically independent configurations), 
(3 = 2.50 (50 configurations), and on 28 4 lattice 
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at 13 = 2.55 (50 configurations), (3 = 2.60 (25 
configurations). To fix the physical scale we use 
the data for the string tension in lattice units 
and use y/a = 440 M eV. 

2. MAIN DEFINITIONS 

To extract P-vortices we use the indirect maxi- 
mal center projection defined as follows 0]. First, 
we fix the maximal Abelian gauge by maximizing 
the functional Rjj(i) = Y^x,n ^[Ux^U+^a 3 ] us- 
ing the simulated annealing algorithm 0. Then 
we replace link matrices U xfl by the Abelian link 

u 11 

variables U xfl -> u xfl = j-^tj. Further gauge 
fixing is done by maximization of the functional 



I? 



Z(2) 



= J2 x ,^ e u x^] 2 - Finally, the maxi- 



mal center projection is performed by replacing 



u xfi -> Z xfl = sign Re u x ^. 



Since we are working on the symmetric lattices 
we have no apriori defined time axis. But when 
we introduce the RxT Wilson loop in the i,4 
plane we define the time direction (/j, = 4) and it 
becomes possible to define the lattice electric and 
magnetic field strength as follows: 



{E? at {x)f = Y,(l-P^)) , 

i=l 



(i) 

(2) 



where plaquette P^ v (x) = }f^dJ x ^ v in case of 
SU(2) fields. The Euclidean action density and 



2 



energy density are defined as follows: 



cr(x) = 

e(x) 



{E? at {x)f- 
-{Et at {x)f 



{Hf at {x)f 
{Hf at {x)f 



(3) 



In the Z{2) projection the static quarks are repre- 
sented by Z(2) Wilson loop, Wz(2), constructed 
from Z(2) links. P^ v {x) entering JQ), (j2J are now 
the average of four Z(2) plaquettes in the {/i, v} 
plain which share the site x. To measure the ex- 
pectation value of the operator 0(x) in the pres- 
ence of static quarks the following quantity should 
be calculated: 

< 0(x) > w = < V M ; > <0(x)>, (4) 

where W is the Wilson loop. In ref.|S] this expres- 
sion was used to calculate the expectation values 
of Z(2) plaquettes in the vicinity of the confining 
string. 



3. NUMERICAL RESULTS 

The results are presented in Figures Q ( a )>(b). 
The Coulomb peaks are not seen for the action 
density. For the energy density one can see rather 
narrow peaks which is in accordance with the re- 
sult of ref. 0] where it was found that the static 
potential in the Z(2) projection is linear down to 
very small distances. 

We define the radius of the Z(2) confin- 
ing string fitting the transverse profile of the 
action density by the expression: c(rj_) = 



Aexp 



(ft)' 



where A and ro are the fitting 



parameters. We treat ro as the radius of the con- 
fining string. To determine the dependence of 
the radius ro on the distance ry between the test 
quark and antiquark we used the data for four dif- 
ferent f3 and various sizes of Wilson loops. Sim- 
ulations at various (3 yield results which are in a 
good agreement (see Fig.0), thus our results cor- 
respond to the continuum limit. To describe the 
dependence of ro on ry we used the fitting func- 
tion ro(ry) = aln(ri|/l fm) + b. The results are: 
a = 0.11(1) fm, b = 0.33(1) fm, x 2 /n dof = 4.31. 
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Figure 1. < a(x) >w (a) and < e(x) >w (b) de- 
fined by eqs. 10, (@J. Calculations are performed 
for /3 = 2.5. The size of the Wilson loop is 8 x 10. 
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Figure 2. The confining string radius as a func- 
tion of In r ii . 
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4. 2q2q system 



The action and energy distribution for four 
static quarks in the SU(2) gluodynamics has been 
studied in 8 . In Fig[3|we show the Z(2) action 
density for two mesons formed by two static qq 
pairs. If the distance between quarks in a meson, 
r* || , is smaller than the separation of two mesons, 
r±, one can see two flux tubes (see FigUJa)). 
When rii = r± there is no any sign of the flux 
tubes (see FigJSJ(b)), the fluctuations are large. 
When r± becomes smaller than rii we observe 
that the flux tubes are formed in the direction 
perpendicular to the initial one (see Fig|3| (c)). 
Thus we observe the effect of the reconstruction 
of the flux tubes. 
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Figure 3. The Z(2) action density for 2q2q system 
at P = 2.5. The size of the Wilson loops is 8 x 10. 



